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The role of relativistic corrections associated with lower 
Dirac components of the deuteron wavefunction, is exam- 
ined for parity-violating (PV) electron scattering. The 
relation between these corrections and negative energy 
components of the struck nucleon in relativistic PWIA is 
elucidated. The model dependence induced by describing 
such effects, using different deuteron vertex functions and 
prescriptions for the half-off-shell nucleon vertices, is com- 
pared against the precision required for studies of nucleon 
strangeness in quasielastic PV e — d scattering. 



I. MOTIVATION 

Measuring the parity-violating (PV) longitudinal 
asymmetry in electron scattering offers the possibil- 
ity to extract the hitherto unknown vector, G S E M , 
and (though to much less extent) axial vector, G S A , 
strangeness nucleon form factors. Several such mea- 
surements will be performed in medium energy fa- 
cilities, such as MIT-Bates [0, Mainz-MAMI || and 
TJNAF g. The first phase of the SAMPLE experi- 
ment Q at MIT-Bates is currently under way and will 
attempt to extract G S M from low momentum trans- 
fer, q 2 = —0.1 (GeV/c) 2 , backward angle, elastic PV 
e — p scattering. However, strong correlation with 
the neutral current (NC), axial, isovector form factor, 
G^ =1 , whose radiative corrections have a large theo- 
retical uncertainty in the Standard Model ||, places 
constraints on the effectiveness of the SAMPLE mea- 
surement. Hadjimichael et al. suggested || that 
a complimentary measurement of the PV asymme- 
try in quasielastic scattering from deuterium, may, in 
conjunction with the proton measurement, allow one 
to determine both G S M and G T A =Y . The rationale of 
this suggestion is that, although both G S M and G^ =1 
appear in the (predominantly transverse) backward 
asymmetry multiplied by magnetic (G^™) form fac- 
tors, they come with different weighting for a nuclear 
target, i.e., (ZG p M +NG n M ) vs. (ZG p M -NG n M ), being, 
respectively, isoscalar and isovector. Thus, for deu- 
terium, this correlation is suppressed relatively to the 
proton by approximately (/ip—/i„)/(/i p +/i„) > 5. Fol- 
lowing this analysis, SAMPLE (in its second phase) 
as well as the E91-017 experiment at TJNAF {§, and 



possibly a future experiment at Mainz M will perform 
such measurements on deuterium. 

In attempting to extract such subtle subnuclear ef- 
fects from a nuclear observable, it is imperative that 
the nuclear physics modeling can be trusted to the 
1-2% level. In Ref. || the theoretical uncertainty was 
estimated for the deuteron case, and was concluded 
that relativistic, rather than rescattering effects, are 
the main source of model dependence in the PV asym- 
metry for quasifree kinematics at |q| > 0.7 GeV/c, 
especially at forward angles (6% effects). This follows 
from a comparison of a traditional potential calcu- 
lation including final-state interactions (FSI), where 
the current operators were truncated to first nontriv- 
ial order in a nonrelativistic expansion, with a plane 
wave impulse approximation (PWIA) calculation, in 
which a nonrelativistic momentum distribution is con- 
voluted with a relativistic off-shell single nucleon ten- 
sor obtained by a generalization of de Forest's "cc" 
prescriptions B to be used with the weak neutral 
current. At low momentum transfer, where relativis- 
tic corrections are -presumably- not important, the 
FSI calculation should be used 0. At TJNAF ener- 
gies, though, the nonrelativistic expansions in the FSI 
calculation are inadequate, while, on the other hand, 
rescattering effects are expected to be suppressed — 
and therefore the PWIA calculation should be pre- 
ferred. In adopting this point of view, the 6% ef- 
fects at |q| > 0.7 GeV/c, which, if taken at face value 
threaten the interpretability of PV deuterium mea- 
surements, are not viewed as a fair measure of nuclear 
uncertainty; rather, the FSI model should be consid- 
ered not applicable for such kinematics 

The PWIA calculation, however, is a semirelativis- 
tic, rather than a covariant calculation, in that it 
treats the deuteron as a nonrelativistic bound state. 
On the other hand, much work has been devoted 
to relativistic treatments of e — d scattering, in the 
sense of consistent nonrelativistic expansions includ- 
ing wave function and boost effects jl(|, as well 
as covariant solutions of the bound state using the 
Bethe-Salpeter equation ]Tl|-[l^] or its "spectator" 
version [ [l4||l5| ] . Given the availability of such relativis- 
tic treatments, and the degree of precision required for 
strangeness studies, the basic objective of this work 
is to calculate the PV e — d asymmetry in a covari- 
ant model. In doing so, our motivation is different 



from most studies of relativistic effects, where observ- 
ables exhibiting large such effects are sought after, in 
the hope to probe details of nuclear dynamics. In our 
case we wish to be as independent as possible from nu- 
clear dynamics. We wish to quantify that the residual 
relativistic effects which are associated with a rela- 
tivistically described deuteron bound state, and are, 
therefore, not included in the semirelativistic PWIA, 
are small and, for the PV asymmetry, less than 1% for 
the kinematics of interest. Since the high-momentum 
bound nucleons are suppressed due to Fermi motion, 
it is reasonable to expect such residual relativistic ef- 
fects to be small. Establishing such a result is im- 
portant for high q 2 quasifree experiments at TJNAF 
with heavier nuclei as targets, for which a semirela- 
tivistic PWIA is possible but a covariant calculation 
is not (notwithstanding the fact that the effects re- 
ported here for the deuteron will underestimate those 
for such heavier nuclei due to the uncharacteristically 
loose binding {Uf ~ 60 MeV/c) of the deuteron). 

Although ultimately we would like to have a covari- 
ant formulation of PV e — d scattering following the 
more fundamental approach of Hummel and Tjon []l2| , 
a simpler approach is taken here, following that of 
Keister and Tjon |ll| and Beck et al. |ll|. It amounts 
to a covariant, or "relativistic" (R) PWIA. Unlike the 
semirelativistic, or factorized (F) PWIA, in which the 
deuteron is treated nonrelativistically, in this case the 
hadronic tensor does not factorize in the strong sense, 
i.e., cannot be written as a convolution of a single- 
nucleon tensor with a momentum distribution. This 
is due to the off-shellness of the propagating nucleon. 
By separating the forward (positive energy) and back- 
ward (negative energy) components of the struck nu- 
cleoli's propagator, we find a "weak" form of factor- 
ization, where the momentum distribution and the 
single-nucleon tensor become density matrices. As a 
byproduct of this construction, we can identify the ad 
hoc prescription in the factorized (F) PWIA of de For- 
est, simply, as a projection of the covariant hadronic 
tensor of the RPWIA onto the positive energy sector, 
which shows the relation between off-shell and rela- 
tivistic effects. Since it is off-shell effects that spoil 
factorization in RPWIA, we expect, and indeed ob- 
serve numerically, that the effects of the projection 
are minimal on the quasielastic peak, where kinemat- 
ics are "maximally on-shell" . 

The structure of this article is as follows. In Sect. II 
we discuss the relativistic calculation (RPWIA) used 
in this work. In Sect. Ill we review the factorized 
PWIA, and discuss the projection by which the FP- 
WIA is obtained from RPWIA. In Sect. IV we show 
numerical results for the cross section and the PV 
asymmetry for the kinematics of interest, and discuss 
the role of the residual relativistic effects in the inter- 
pretation of experiments designed to study nucleon 



strangeness. Our conclusions appear in Sect. V and 
some details of the calculation in the Appendices. 



II. THE HADRONIC TENSOR IN 
RELATIVISTIC PWIA 

The kinematics for deuteron electrodisintegration, 
in PWIA and in the target rest frame, d M = (Md, 0), 
are shown in Fig. 0(a). The momentum transfer is 

= (lu, q), and, for the moment, we assume that the 
nucleon p^ = (E p , p) is detected in coincidence with 
the electron e' M = (e', u'), while the other nucleon 
= (E n ,—k) is not. Here d 2 — Mj, p 2 = n 2 = m 2 . 
The cross section is obtained from 

(h) _ a? du'dp dn rn* m 2 4 

~ q+ (27r) 3 Md ee , EpEn d W W^U , 

(1) 

where = d^ + q^, P^ = p^ + n^, is the leptonic, 
l»v = Y] < e', h'lj^e, h >*< e', h'\^ v \e, h > , (2) 



and the hadronic tensor 



J2 . <f\Jg\i>*<f\JZ\i> 

' ^snins c c 



(3) 



and where spin d.o.f. have been averaged over in the 
unpolarized target (\i >= \d >) state and summed 
over in the final electron and np states. The plane- 
wave impulse approximation (PWIA) amounts to (I) 
using one-body current operators J M , (II) disregard- 
ing final-state interactions, and (III) assuming that 
the detected nucleon is identified with the (off-shell, 
k 2 < m 2 ) nucleon k^ — (Ek , k) struck by the virtual 
boson P,|l7|. Under these assumptions the PWIA 
amplitude reduces to that of Fig. |l|(b). The inclu- 
sive cross section in PWIA with polarized electron 
beam and unpolarized target is obtained by integrat- 
ing over G? 3 n and dfl^ using the energy-momentum 
(5-function JOI 

dV 1 ) 1 1 f Y m 2 

dn^ = ^M' d ( ™ 4, d|k| |k| 9 ^ 2 + |k| 2 

x (v L W l L + v t Wt + h v T ,W^,) (4) 



= a M \ v L R + v T R + h v T 'R 



T' 



Here om is the Mott cross section. The limits of in- 
tegration, \y\ and Y, as well as the responses R are 
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functions of |q| and u. 
included in the factors 



The lepton kinematics are 
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The incoherent sum over nucleons in Eq. ( |) is a con- 
sequence of assumption (III) above [0 . The respones 
W l are obtained as components of the hadronic tensor 
in Eq. (§. 





(a) (b) 
FIG. 1. Deuteron electrodisintegration in PWIA. 



In relativistic (R) PWIA the hadronic tensor is ob- 
tained from the amplitude in Fig. |l|( 



< f\j»\i > ~u(p)r»-t^X v Cv(n) . (6) 

Here is the half-off-shell nucleon vertex, £„ is the 
deuteron polarization vector and X v is the d-p-fi ver- 
tex. Having applied charge conjugation to the on-shcll 
nucleon \n >, rather than the off-shell one \k >, X p is 
obtained from the Blankenbecler and Cook ver- 
tex, V M , according to 



x v = C T VjC T 
B 

= A-y v H r v 

m 



2m 



G 

F lv + —r u 
m 



(J) 



with C the charge conjugation operator, and r v = 
(n v - k v ) /2. The vertex functions A, B, F, G depend 
on k 2 and can be written in terms of the 3 S\ (w), 
3 Z?i (w), 3 Pi (vt) and 1 P\ (v s ) relativistic deuteron 
wavefunctions [jll|l4|. In the on-shell matrix element 
u(k)X„v(n) the off-shell F, G vertex functions would 
not contribute. Squaring the amplitude in Eq. (||), 



1 the T' response is included anticipating the discus- 
sion of parity-violation below; it is identically zero for 
purely electromagnetic (EM), inclusive, unpolarized tar- 
get scattering. 



summing over spins in the final state, and averag- 
ing over deuteron polarizations, one obtains for the 
hadronic tensor in RPWIA 



RPWIA 



5pA + d p dx/Mj \i 



3(fc 2 -m 2 ) 2 



Tr 



2m 



X p (Ijk+m) 



2m 



(8) 



where X v = 70^70 and similarly for f M . Following 
Keister and Tjon [Oil we define auxiliary vectors "P M , 
and scalars f3, T through 

= (- 9p x + d p d x /Ml) X x (i-m)XP (9) 



2m 



2m 



(10) 



implying 



(2m) 2 S^ = 2 [m/3 + (k-V)} k» + (m 2 - k 2 )V^ 
(2m) 2 T = (k 2 + m 2 )(3 + 2m(k-V) . (11) 



We thus obtain 



4m 

"RPWIA - Q fi,2 ^2\2 



3(k 2 — m 2 



Tr i £±I f> t+IH r* 



2m 



2m 



(12) 



Expressions for V, (3, and the hadronic tensor can be 
found in Appendix A. Equations ( p"2| ) and ( |TT| ) provide 
a reference point for a comparison with the factorized 
(F), or "semirelativistic" , PWIA, that is widely used 
in intermediate energy physics. 



III. REDUCTION OF RPWIA TO FPWIA 

A. Factorized PWIA revisited 

In FPWIA the cross section is written as a con- 
volution of a nonrelativistic momentum distribution 
p = /o(|k|) with a single nucleon cross section [S[l8) 



rf 2 qW 
dfl P dw 
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2n (7m 



d\k\ |k| P (|k|) — 

y\ QV 171 
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k| 2 

hv T >W T ,\ , (13) 



The single nucleon responses Wl,t,t' are obtained 
as components of a covariant single nucleon tensor 
Tt^ , describing the interaction of a moving, off-shell 
nucleon with the virtual photon. The most general 
form of the off-shell vertex has been obtained by Bin- 
cer [^lj . The special case of the half-off-shell vertex 
u(p)r M involves four independent vertex functions of 
two scalar variables, e.g., q 2 and k 2 , and collapses 



to the familiar form involving just two form factors 
(Dirac, F\, and Pauli, F 2 ) when operating on an 
on-shell spinor u(k). The half-off-shell vertex func- 
tions of composite particles cannot be uniquely de- 
fined |^2|,^3| . They can of course be computed in spe- 
cific models 24 2q|, but then the whole reaction (i.e., 
the self energy and the d-p-h vertex) should be, as 
well. Since such simple models are not realistic (i.e., 
they do not reproduce the observable on-shell form 
factors F\ and F2) it has become customary to resort 
to ad hoc prescriptions for treating the half-off-shcll 
single nucleon tensor, the most widely used of which 
are the "cc" prescriptions of de Forest |gj . These pre- 
scriptions involve a set of operations regarding 

I. the wavefunction of the struck off-shell nucleon. 

II. the form of the half-off-shell vertex. 

III. the treatment of current conservation. 



With respect to III, one notices that, in the half-off- 
shell case, the one-body current is not conserved 
since q IJ ,u(p)j tl — u(p)(m — The de Forest "cc" 
precriptions enforce current conservation J • q = Jqlu 
by eliminating the third component of the current in 
favor of the charge density, where q defines the z-axis. 
Alternatively, one may choose to eliminate the charge 
in favor of J3, or not enforce current conservation at 
all. In Ref. [^6| it was shown that the treatment of 
current conservation induces, especially in high miss- 
ing momentum kinematics in (e, e'p) processes, a sub- 
stantial uncertainty. Recently, the three ways of treat- 
ing current conservation have been shown to be equiv- 
alent to specific gauge choices within a family of co- 
variant gauges |27j| . Thus, the RPWIA approach de- 
scribed here is covariant, but not gauge invariant. 



B. Strong vs. Weak Factorization 



With respect to (I), the struck, off-shell nucleon 
wavefuction is taken to be that of the correspond- 
ing "would-be-on-shell" spinor u(k), with same spa- 
tial momentum k, but (positive) energy E% = E n = 
\Jm 2 + k 2 > Ek- As we discuss in the next section, 
by reformulating (I) in a more transperent fashion, it 
is this approximation that differentiates between RP- 
WIA and FPWIA. 

With respect to (II) the following two vertices are 
employed 



1 2 



F17" 



i^a^F 2 
2m 



(f 1 + W-^+7^^ 
2m 2m 

(Fx+F 2 )^-^±^F 2 , (14) 
2m 



where the second equation in (|l4|) is to be understood 
in the sense of the Gordon decomposition of u^T^. 
In the vertex appears the auxilliary four-vector 
k^ = (E k ,k), instead of = (E k ,k). In the on- 
shell case, the Gordon decomposition of u(p)T2u(k) 
is complete and the T^, vertices matrix-element- 
equivalent. The difference between results obtained 
using the two forms of the vertex is interpreted as 
an estimate of the theoretical uncertainty induced in 
the modeling of the reaction because of our inability 
to calculate reliably the half-off-shell vertex functions 
from first principles flpq]- Thus, under I and II, the 
single nucleon tensor reads 



"fpwia — I 2m 1 ' 2 2m 



1 1,: 



(15) 



The question we address here is under which set 
of operations (related to assumption (I) above) is the 
FPWIA of Eqs. @ and |jj) obtained from the RP- 
WIA of Eqs. (|) and (ITT 



#rpwia ^ (^) 3 M d P Kpwia • (16) 
We start by observing that from Eq. (nw one has 



$ + T = 2[mf3 + (k.V)]^ + ?^[T?-f3} . 



Am 2 



(17) 



Since V v is obtained from the deuteron vertex func- 
tions [cf. Eq. (||)], the second term in Eq. (jl^) does 
not allow to factor all dependence on details of the 
deuteron break-up amplitude out of the trace of Eq. 
(|l2|), thus defining a "single nucleon" tensor. How- 
ever, if one ignores the off-shellness of the struck nu- 
cleon and sets k^ — > fc M in Eq. (pTl), one obtains 



(18) 



4to 2 (#+T)^2 m/3+(k-V) fJS+m) 



Comparison with Eqs. (gj), @, Q and @ shows 
that the cross section has indeed the factorized form 
of the "semirelativistic" FPWIA, provided that the 
nonrelativistic momentum distribution p is related 
to the quantity [to/3 + (k-V)] / [l2n 3 M d (k 2 - m 2 ) 2 ] . 
We shall return to this later on. The relevant obser- 
vation at this point is that the RPWIA is not factor- 
izable, unless the off-shcllncss of the struck nucleon is 
ignored. To probe this issue further, we use the fa- 
miliar technique of interpreting the propagation of the 



where (I) has resulted in the positive energy projector 
with the auxilliary four- vector k. 



analogous observations can be made for the case of deep 
inelastic scattering processes p8(. 



A 



off-shell struck nucleoli as a superposition of positive 
and negative (on-shell) energy states (see, e.g., Ref. 
p9[ , or Arnold et al. p5[ , where a similar formalism 
is developed for elastic e — d scattering) 



1/tA-m 
2m 



2m 



/- 



— m 
2m 



where 

/± = (E k ± E k )/2E k , and fc£ = (E k ,±k) 



(19) 



(20) 



The decomposition, Eq. (|19|), is frame-dependent. 
Notice, in particular, that in the deuteron rest frame 
the on-shell states are the struck and spectator nu- 
cleon, respectively, i.e., k^_ = k^- and kt = We 
now write Eq. (H9) in terms of the direct products of 
these on-shell states 



ftA-m 
2m 



= /+ E M ( k ' s )""( k ' s ) + f~ E v (~ k > s M~ k > s ) 

S S 

(21) 



and perform this decomposition on Eq. (|8|). In the 
shorthand notation x+ = u (k, s) and \ a _ = v(— k, s), 
the trace of Eq. (|lj) now reads (Dirac indices explic- 
itly shown and implicitly summed over) 

~E E fah\V+P]^Xl]^a] k [^] kl \HM lm 
s,r a,h=± 

xrutxants]* ■ ( 22 ) 

By inspection, this can be cast in the form 

^.PWIA = (^) 3 Md E E P<*T. Kirs • (23) 
spins a,b— ± 



(A; ± -"P)±m/3 



P+-;rs = /+/- 

= P- + ;rs 



127T 3 Af d (fc 2 - TO 2 ) 2 

Ek{t.,r'V)-[{k+-P)+m0\ 



(27) 



l2ir 3 M d {k 2 - m 2 f 



(28) 



Thus, the diagonal terms (+, +) and (— , — ) are diag- 
onal in spin indices as well. Setting 

P±±;rs = P±± Sr,s, E §r ' s W ±±;ra = ^±± i ( 2 9) 
spins 



we obtain for the diagonal contributons 
P±±-rs n ±± . rs — p±± n ±± , 



(30) 



where we make use of J2 s £ss = to identify 

[cf. Eq. @], P++ = f*[mp+ (k-V)} 



rl FPWIA 



/l2n 3 M d (k 2 — m 2 ) 2 , and analogously for the 



) 



contribution. The observation made in Eq. (JL8J) can 
thus be formally stated as follows: the first step (ac- 
tually the only step, as shown below) required for the 
reduction to FPWIA [cf. Eq. ([l6|)] is a projection to 
(+, +) energy sector. With respect to the off-diagonal 
(in energy indices, a ^ b) contribution, notice that the 
relevant momentum distribution does not appear to 
be diagonal in spin space either. The second term 
in Eq. (|l7|), which, as noted above, spoils strong fac- 
torization, is an off-diagonal energy (+,— ) term, as 
suggested by writing 



rn 



k 2 = E 2 ~ E 2 = -2/+/_ [m 2 + (k-V)} • (31) 



This is verified by explicit evaluation in Appendix B. 



A weak form of factorization, therefore, still holds, 
with momentum distribution p a b,rs now being a den- 
sity matrix in both spin and energy projection indices 



Pab;rs I a lb 



m 



6ir 3 M d (k 2 - m 2 ) 



(24) 



C. Nonrelativistic Momentum Distribution 

We now address the relationship between p ++ and 
the familiar, nonrelativistic deuteron momentum dis- 
tribution p(\k\) ~ [w 2 (|k|) + w 2 (|k|)] . In terms of the 
wavefunctions u, w, v s and v t , the vertex functions in 
Eq. (0) are given by [pi 



and similarly for the single-nucleon tensor 



2m 



Using v{— k, s) = 7 5 7°u(k, s) and (30 1 



u(k, s) u(k, r) 



ft S s ,r +7 5 



2m 



where £g r = u(k, r)7 5 7 M it(k, s), we find 



(25) 



(26) 



A = Af 



G 



B = Af 



m m(2E k + m) 
-u A = w ■ 



E k + m" ■ V2|k| 2 

2Mm 2 [ u E k + 2m 

M d [E k Arm ~ V2|k| 2 

3 m 



- w 



F = -2^ v / 2Md^y| 



Ikl V 2 Vt 



2Try/6M~ d 1 ^ T v s 
k 



(32) 



where AT — Ti\/2Md{E k — -Efe). Intuitively, we expect 
that p ++ reduces to p in the nonrelativistic limit, that 
is, by ignoring the P states, and, additionally, making 
a (|k|/m) 2 expansion |n]. It turns out, however, that 
neither approximation is necessary, and, in fact, p ++ 
and p are identical. To see this, first notice that the 
overall contribution of the P states must cancel out 
once the projection to the (+, +) sector is made, since 
such states represent "lower" Dirac components ip~ of 
the deuteron wavefunction: from Eq. (38) in Ref. 14 
we have (in our notation) 



3 5i, 3 Di G tl>+~x+{X-Ov(n) 
"Pi, 1 Pi G ^-~x-(*-0«(n) 



(33) 



By construction, the (+, +) sector does not contain 
X- terms. Thus, p ++ is oblivious to the P states 
(we have verified this numerically, as well), which can 
therefore be neglected for the purposes of writing p++ 
in terms of wavefunctions. Notice that, although in 
this case F — > 0, there is still an apparently nonvan- 
ishing contribution to p ++ from the "off-shell" vertex 
function G; however, this can be show n to v anis h iden- 
tically. For example, from Eqs. (|A2| ) and ( |A3|) in the 
appendices, and using Eq. (32), we have for the G 2 
terms in the deuteron rest frame 



0{G 2 ) ~ m 2 + k 2 + (k 2 -m 2 ) 



(k-n) 



(k-k) 



- 2{k-k) 



+2(k-n) - 2(k-n) 
+ E 2 + 2E k E k + (E 2 -E 2 k 



E 



2 -'k| 2 



-jL(E 2 E 2 -\k\ 4 )-2E k E k + 3\k\ 



TO 2_ 2 ^| + ^2 +3|k| 2 +2 ^r 



= 



(34) 



and similarly for the AG and BG terms. Thus, once 
the P states have been shown to decouple, the nonva- 
nishing contribution comes f rom the A and B vertex 
functions alone. From Eq. (|A2| ) and Eq. (A3) with 
V 2 = -Ikl 2 , we find 



m 



(3 + k-V = A 2 (3m 2 + 3(n-k) - 2(N-k)) 



V 2 

+ B 2 — (m 2 - (n-k)) - 2AB ((V-n) - (V- 



+I2m 2 [u - w/V2j w/V2 
= 67V 2 m 2 (u 2 + w 2 ) . (36) 

After some suppressed algebra we find that retain- 
ing all terms modifies Eq. ( |36| ) multiplicatively, by 
E 2 /m 2 . Thus, we finally obtain 

m/3 + k-V = 6M 2 El (u 2 + w 2 ) 



=> p++ = n 



mf3+{k-V) 
l2ir 3 M d (k 2 - m 2 ) 2 



6M{E k + E k ) 2 E 2 {u 2 +w 2 ) 
48ir 3 M d E 2 (k 2 -m 2 ) 2 



-in 



where we have used Eq. (20) and k 



(37) 



El 



A 2 (6m 2 +4|k|- 



-AAB\k\ 2 + 2B 2 Hr 



p++ is precisely the nonrelativistic momentum dis- 
tribution, apart from an overall normalization, since, 
when P states are present, p++ is not suitably nor- 
malized, i.e., J °° dKK 2 [u 2 + w 2 ] < 1. Before proceed- 
ing to the next section, two comments seem appro- 
priate. First, it is important that p ++ contains k-V, 
rather than k-V. This can be seen by observing that 
the leading order coefficient of the B 2 term must be 
~_|k| 4 , otherwise p ++ diverges as |k| — > [cf. Eq. 
©]. Using k-V would add to the RHS of Eq. © a 
term B 2 \k\ 2 E k (E k — E k )/m 2 , which, because of the 
deuteron binding energy b, would diverge as |k| — > 0, 
since E k - E k = \b\ + 0(|k| 2 /m 2 ). Second, we em- 
phasize that Eqs. (^J) and ( J35| ) are only valid after 
the overall contribution of P states has cancelled out. 
In other words, the F and G vertex functions do con- 
tribute in p ++ when the P states are present; their 
contribution, though, serves to cancel that of the P 
states from the A and B terms. 



IV. VALIDITY OF THE REDUCTION 

A. Kinematical considerations 

In the previous section it was shown that the de 
Forest type of factorized PWIA is obtained from the 
RPWIA by projecting onto the the positive energy 
sector of the propagating struck nucleoli. The legiti- 
macy of this projection depends on / + ~ 1, /_ << 1: 
the more off-shell is the struck nucleon, the less valid 
the projection. Energy/momentum conservation at 
(3Jt)ie d-p-n vertex in the deuteron rest frame implies 
that the struck nucleon is off-shell by 



To leading nonrelativistic order we find from Eq. ( ]35| ) 
and Eq. (132) (without u t !) 



E k -E k = 2V™ 2 + k 2 -M d 



m/3 + k-V =Af 6m 2 [u-w/v2 



< dm 2 w 2 



k 2 

\b\ + ~ 

m 



(38) 



r; 



where the ellipses denote higher order relativistic 
terms. The binding effect is negligible for the 
deuteron « 2.22 McV) and |k| is effectively 
restricted by the small deuteron Fermi momentum 
(kp ~ 60 MeV/c), resulting in a rapid falloff of the 
momentum distribution p(|k|). The projection seems, 
therefore, meaningful (at least in the deuteron rest 
frame). Some caution must be exercised, however, for 
two reasons. First, there may be cases, where, despite 
~ 1, backward propagation is essential. An exam- 
ple is low energy Compton scattering from, e.g., a 
proton, where, to leading order, forward propagation 
cancels between the s and u channels and the classi- 
cal (Thomson) limit is entirely due to the backward 
propagation terms |29). Second, the applicability of 
the Fermi motion argument depends on the kinemat- 
ics (choice of energy transfer to for fixed momentum 
transfer |q|), which determines the integration region 
for inclusive scattering [cf. Eq. ((u|)]. Overall en- 
ergy/momentum conservation implies |lq,BQ| for the 



integration limits k„ 



\y\, K 



Y 



1 m 2 



y (\q\, w ) = (M d +w)^-- — -if, (39) 

and Y = y + |q|, with s = (M d + ui) 2 - |q| 2 . On the 
quasielastic peak (u> — u> qe ) 



M d + uj qe = yjm 2 + |q| 2 4 
2m 



J qe. 



(40) 



the lower limit \y\ = 0, and, therefore, the contribu- 
tion of small momentum components |k| saturates the 
nuclear response in Eq. ( |T3| ) . Higher momentum com- 
ponents contribute as well, however, their importance 
(as weighted by the rapidly falling momentum dis- 
tribution), is negligible. Away from the quasielastic 
peak, though, y ^ 0. As an extreme example consider 
the situation at threshold, where s — Am 2 , and thus 



Md + U t hr 



yjAm 2 + |q| 2 

j^+0(\b\) 
4m 



(41) 



In this case \y\ = |q|/2, and the struck nucleon will be 
off-shell by 1 - E k /E k > (2 - M d /y/m 2 + y 2 ), which, 
near threshold, translates into 24%, 63% and 115% 
for typical TJNAF processes with |q|=l, |q|=2 and 
|q|=4 GeV/c, respectively. 



B. Electromagnetic responses and the cross 
section 

The effects of the projection on the longitudinal 
(Rl) and transverse (Rt) responses entering the in- 
clusive cross section can be seen in Figs. ^-|| The 



responses are decomposed into (+,+), (+, — ) and 
(— , — ) contributions, after carrying out the spin sum- 
mation in Eq. (p3|). By (+, — ) we denote the total off- 
diagonal contribution, i.e., (+, — ) plus (— ,+) in the 
notation of Eq. (^). Unlike the (+,+) and (-,-) 
contributions which are always positive in the case of 
the electromagnetic responses, the off-diagonal ones 
appear with either sign. For numerical evaluation, it 
is easier to project directly from Eq. ( [12] ) by apply- 
ing the decomposition (19) to Eq. (|l0|), and then use 
f + — /_ = 1 to arrange the resulting expression in 
/+) /+/- an d f- terms. We have verified that this 
is identical to the direct trace evaluation of Eq. fl23| ) 
described in Sect. Ill and Appendix B. 

The residual relativistic effects we are interested 
in correspond to the difference between RPWIA and 
FPWIA results. From the discussion of the previous 
section follows that the R ++ responses are identical 
with the FPWIA ones modulo an overall normaliza- 
tion, since J cIkk 2 (u 2 + w 2 ) < 1 when P states are 
present, while in FPWIA one uses wavefunctions nor- 
malized to unity. This results to a maximum 1.4% 
effect (with the A = 1 vertex functions, see below). In 
the case of the PV asymmetry, where such 1% differ- 
ences are important, this renormalization is irrelevant, 
as it cancels in the ratio of cross sections. Thus, in 
the following we use the terms (+, +)-projected and 
FPWIA interchangeably. 

We present results for |q| = 0.3 GeV/c, represen- 
tative of a low energy experiment such as the SAM- 
PLE project, and |q| = 1 GeV/c, representative of a 
medium energy experiment at TJNAF. Current con- 
servation has been enforced ( "cc" ) by eliminating J 3 , 
as discussed in Sec. III. For the on-shell form factors 
appearing in the electromagnetic vertices 2 [cf- Eq. 
(|l4|)], we assume that they will be experimentally well 
known by the time the PV measurements will have 
been completed, and thus use the parametrization in 
Refs. |3^,|33| rather than model calculations. In this 
parametrization G P E = (1 + 4.97t)~ 2 , G^j" = fi P:n G p E 
and G% = -rG' A V(l + 5.6t)-\ with t = \q 2 \/Am 2 
and where fj, p — 2.79, /!„ = —1.91 are the proton and 
neutron anomalous magnetic moments, respectively. 
For the four vertex functions in the d-p-n vertex we 
use the spectator equation solutions obtained by Buck 
and Gross 10. These solutions are labelled by a 
parameter A that interpolates between pseudovector 
(A = 0) and pseudoscalar (A = 1) 7riVW coupling 



A-z(l-A)-^ 
2m 



(42) 



From the interpretation of the P states as lower Dirac 
components of the deuteron wavefunction, one ex- 
pects that pseudocsalar (ps) coupling would enhance 
the role of the P states, since 75, being off-diagonal, 
mixes lower and upper components. As discussed in 
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Ref. H] this is indeed the case, and the importance 
of P states increases rapidly between A — (pv) and 
A = 1 (ps). 




FIG. 2. Percentage decomposition of inclusive elec- 
tromagnetic responses (left panels, Rl, right panels, Rt) 
using the Fj vertex, at |q| = 300 MeV/c. Dashes: 
(+, — ) contribution; dot-dashes: ( — , — ) contribution; 
solid: total contribution of negative energy components, 
i.e., 1 — R ++ /R. Results shown for different choices of 
deuteron vertex functions. Fhe solid triangles mark the 
quasifree region [cf. Eq. (0)]. 

Commensurate, looking at the "BG" panels in Figs, 
we see that the difference between FPWIA and 
RPWIA, i.e., the error induced by projecting onto 
the (+, +) sector, significantly increases with A. It 
is well known, however, that the ps and pv irNN La- 
grangians are related via a field redefinition (to lowest 
order in A, at higher orders the field transformation 
generates contact terms, as well). According to basic 
field theory properties, the amplitude of the overall 
on-shell process should not depend on A f|4|,|23| ■ The 
fact that it actually does, partially reflects the incon- 
sistency of computing the d-p-ii vertex functions in 
a specific model (both the on-shell, A, B and off- 
shell, F, G) ones, while not using the same model 
for the vertex functions of the photon-nucleon ver- 
tex. In particular, it is inconsistent to use the off- 
shell vertex functions F and G in the d-p-ii vertex of 
Eq. (|7|), while ignoring altogether the two extra ver- 
tex functions entering the half-off-shell jNN vertex. 



As a side remark, one notices that, while realistic (in 
the sense that deuteron properties can be reproduced 
using them) model solutions have been obtained for 
F and G, first-principles calculations of the off-shell 
electromagnetic vertex functions are indicative rather, 
than realistic Conceivably, this reflects the fact 
that in the deuteron off-shell vertex one needs the 
constituent off-shell, while, in the jNN off-shell ver- 
tex, one needs the bound state itself off-shell, which 
is, apparently, more difficult to describe. 
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FIG. 3. As in Fig. H, but with the vertex 



It appears, therefore, reasonable to ignore the off- 
shell vertex functions in the d-p-ii vertex, as long as 
only on-shell electromagnetic form factors are used. 
This point of view has been adopted, e.g., in Ref. fllrjf , 
where a parametrization of the deuteron vertex by 
Locher and Svarc fjSfl using only the A and B ver- 
tex functions was employed. Since our primary mo- 
tivation in this work is to investigate as widely as 
possible the extent to which effects associated with 
the relativistic aspects of the deuteron vertex modify 
the factorized PWIA, we will, nevertheless, use the 
off-shell deuteron vertex functions. In an attempt to 
partially alleviate the inconsistency of this scheme, we 
also show some results, whereupon we still use only 
two vertex functions for the ^NN vertex, but modify 
the on-shell form factors according to 



s 



F(q 2 ,k 2 ) = F Galater (q 2 



F ps (q 2 ,k 2 ) 



F ps (q 2 ,m 2 ) 



(43) 



< V2fc F |q|/2v^ 



(44) 



allowing for a dependence on the other scalar variable, 
k 2 . For the modification factor we use the one-loop 
calculation with pseudoscalar (ps) irNN coupling of 
Naus and Koch ]2^ ]. Admittedly ad hoc, this pro- 
cedure (first adopted in Ref. |3(|) incorporates, nev- 
ertheless, what has been shown to be - within a toy 
model - the prominent feature of the EM half-off-shcll 
vertex. For consistency, we only apply this modifica- 
tion scheme with the pseudoscalar (A = 1) deuteron 
vertex functions. For comparison purposes, we also 
show some results in the Locher-Svarc scheme (la- 
belled "LS"). We find, in general, that the effects of 
the projection are large in this scheme, comparable in 
magnitude to those with the A = 1 Buck-Gross vertex 
functions. 
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FIG. 4. As in Fig. @, but at |q| = 1 GeV/c. 

Looking at Figs. || and ||, one sees that at the low 
momentum transfer value (|q| — 300 MeV/c), the ef- 
fects of neglecting negative (— ) components in the 
longitudinal response, Rl, are negligible for all en- 
ergy transfer values w considered, independently of 
the deuteron vertices and form of the EM current op- 
erator used. For the transverse response, Rt, more 
significant effects are seen, particularly with the 
vertex. Specifically, with the vertex (Fig. ||) the 
effects increase with A but are still below 2% inside 
the quasielastic ridge 



For the vertex, however, 5-10% effects are seen 
(Fig. ||), especially with the A = 1 vertex functions. 
In all cases, the effects are minimal for quasifree kine- 
matics, as anticipated from the discussion in the pre- 
vious section. At this momentum transfer value, we 
find that a modification according to Eq. d43| ) does not 
produce an observable effect in the the A = 1 panels 
of Fig. | 
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FIG. 5. As in Fig. |§ but with the vertex. In the 
X = 1 panels results are shown with and without the mod- 
ification of the on-shell form factors according to Eq. (kq). 



Unlike the case at the low momentum trasnfer, at 
q = 1 GeV/c (Figs. | and |), R L is affected by the 
projection more than Rt- Inside the quasifree region, 
though, the effects do not exceed 3% for either Rt or 
Rt- Thus, although the behavior of Rt is qualita- 
tively similar to what happens at |q| = 300 MeV/c, 
the effects - particularly with ccl - are, within the re- 
spective quasifree regions, smaller than at |q| = 300 
MeV/c. This behavior does not contradict the re- 
marks following Eq. (|4l]), because threshold is close 
to the quasifree region at |q| = 300 MeV/c (where 
w thr = 24 MeV), but not at |q| = 1000 MeV/c, 
where w t hr = 250 MeV; that large effects are expected 
towards threshold can be seen by extrapolating the 
A = 1 ccl Rt panel in Fig. ||. This behavior in the 
quasifree region verifies that the residual relativistic 



Q 



effects we describe are rather independent of the mo- 
mentum transfer q 2 . If anything, they appear more 
important for low q 2 values, as can be also seen in 
Fig. ^, where the ratio of the (+, +)-projected (i.e., 
FPWIA) to the full (+,+) ffi (+,-) 8 (-,-) (i.e., 
RPWIA) responses is shown, for quasifree kinemat- 
ics. The A = 1 vertex functions are used, in order 
to display the maximum effect. The Rt response ob- 
tained with the EM vertex shows > 2% effects 
below |q| = 600 MeV/c. At the SAMPLE kinematics 
(|q| = 300 MeV/c) the effect on Rf 1 is 5%. 
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(45) 
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where we have set A = k — k = (Ek — Ek , 0) . The 
A terms do not induce a large effect. In fact, they 
do not contribute anything when /i = 0. In the FP- 
WIA, the — m) term does not contribute because 
of the (p+m) projector onto the (+, +) sector. Thus, 
the ccl and cc2 longitudinal responses (both R L and 
its parity- violating counterpart, R\ v ) are identical 
in FPWIA |^,^6| while the transverse ones differ, 
though not by much. This is not discernable from 
Figs. |]-|], but can be seen from Fig. |?], especially 
the inner plot, where the FPWIA cross sections differ 
insignificantly between ccl and cc2 throughout the ex- 
citation energy region (solid line). The (Jji— m) term, 
however, contributes to the (+,— ) and (— ,— ) parts 
of the RPWIA. It is therefore not surprising that a 
significant difference is seen in Fig. between the ccl 
and cc2 cross sections in RPWIA (dot-dashed line in 
inner plot) for non-quasifree kinematics, to be con- 
trasted with the coresponding behavior after (+,+) 
projection (FPWIA, solid line). 
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FIG. 6. The ratio R ++ /R of electromagnetic (L,T) 
and parity-violating (Lav ,Tav ,T va ) responses accessi- 
ble in inclusive d(e, e') scattering, for quasifree kinematics 
(lj — ujqe) as a function of momentum transfer |q|. The 
A = 1 Buck-Gross deuteron vertex functions are used. Re- 
sults are displayed for both ccl and cc2 prescriptions for 
the single- nucleon tensor (the T and Tav curves cannot 
be disentangled in the ccl panel). 
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FIG. 7. Inclusive d(e,e') cross section using the 
covariant (Ri) and (+, +)-projected (Fi) responses, 
as a function of the final state excitation energy 

The A = 1 vertex func- 



E, 



y/(M d + Wf 



2m. 



tions and both i = 1, 2 EM vertex prescriptions are used. 
Ratios between different calculations are displayed in the 
inner graph. Data from Ref. JS^]. 



Returning to Fig. [|, note that the (— , +) terms are 
not always dominant over the (— ,— ) ones (see, e.g., 
the A = 1 panels), as would have been expected given 
the extra suppression of the latter by /_// + . Finally 
the A = 1 results in Fig. [| show an observable, but 
not dramatic change when the EM form factors are 
modified according to Eq. (p3[). 

Turning to a direct comparison between ccl and cc2 
results, we write, using Eq. (14), 



Notice also, from the a^/a cc \ (labelled F\/Rx) 
and <j^2 /°cc2 (labelled F2/R2) ratios in the inner 
graph of Fig. ^, that the negative energy components 
play a larger role when the T 1 ^ vertex is used. In fact, 
this is a common and rather unexpectred feature of 
all Figs. the effect of the projection is more sig- 

nificant when ccl is employed. Given that either EM 
vertex amounts to an equally ad hoc choice, this ob- 
servation cannot justify the overwhelmingly dominant 
use of the 'ccl' prescription in data analyses, if not 
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suggesting that relativistic effects are more appropri- 
ately incorporated when using cc2. 



hadronic tensor corresponding to the A 1 — Az inter- 
ference, 



C. Parity-violating responses and asymmetry 
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(48) 



In inclusive (e, e') scattering of longitudinally polar- 
ized electrons (helicity h — ±1) from an unpolarized 
target, the asymmetry A — (da + — da~)/ (da + +da~) 
is purely parity-violating. There are two sources of 
parity-violation, electroweak (7 — Z°) interference, 
and nuclear parity-violation. In the latter case, the 
weak interaction leads to a parity-violating compo- 
nent in the NN Hamiltonian, which allows parity mix- 
ture with P states (not to be confussed with the rel- 
ativistic P states discussed so far, which correspond 
to lower components of the coupled Dirac wavefunc- 
tions and therefore have the same overall parity as the 
S and D upper components [Q; these P states, on 
the other hand, arise in a nonrelativistic framework 
and have genuinely opposite parity). Nuclear parity- 
violation appears suppressed relatively to 7 — Z° in- 
terference by a/a s , where a s denotes the strength 
of the (parity-conserving) strong interaction. In a 
recent calculation |38| the nuclear PV asymmetry 
appears, indeed, negligible (suppressed by three or- 
ders of magnitude) compared to the asymmetry due 
to electroweak interference over the whole range of 
momentum transfer values considered (q € [0, 1000] 
MeV/c). Thus, in the following we focus on the PV 
asymmetry generated from the interference of electro- 
magnetic (photon-exchange, Ay) and neutral current 
(Z°-exchange, Az) amplitudes. In the one-boson ap- 
proximation, and for momentum transfer values much 
smaller than the electroweak scale, \q 2 \ « Myy, the 
asymmetry is given by |33[] 



A = 



G\q 2 \ a A [vlR^v + v t Rav] + a W T >RvA 



2?rav / 2 



vlRl + vtRt 



(46) 



Here G is the Fermi constant and ay (cm) are the 
vector(axial-vector) couplings in the e — Z° vertex 
(oy7/i + 0^757^). To understand the V — A structure 
of Eq. ([l6|), it is useful to observe that the helicity 
[h = ±1) dependent part of the leptonic tensor 
corresponding to the A 1 — Az interference, reads 



into H^ v = S^ v + iA^ v . For unpolarized targets, an 
antisymmetric part A^ v may only arise from axial- 
vector coupling. From general considerations in in- 
clusive electron scaterring ]l9| ] with polarized beam 
and unpolarized target, the L and T responses are 
generated by contracting the symmetric parts of the 
leptonic and hadronic tensors. Thus, from Eq. ( fl^ ) 
the L and T parity-violating responses in the numer- 
ator of Eq. (^) arise from interference of the lep- 
tonic axial-vector with hadronic vector pieces (hence 
the AV index). The T" response, on the other hand, 
arises from contracting the antisymmetric tensors; it 
does not, therefore, contribute in unpolarized scatter- 
ing hence, docs not enter the helicity sum in the 
denominator of Eq. (^6|) - and has a leptonic/ vector 
with hadronic/ axial- vector (VA) interference charac- 
ter. Thus, the nucleon axial-vector form factor G^ =1 
that plays an important role in motivating these stud- 
ies, as described in Sect. I, appears in the RyA re ~ 
sponse. Notice, though, that since in the Standard 
Model ay = -1 + 4 sin 2 6 W ~ -0.092, while a A = -1, 
the RyA response is effectively suppressed with re- 
spect to the other responses in electron scattering (but 
not in neutrino scattering, where ay = 1 (33]). 

The PV responses are computed from an interfer- 
ence hadronic tensor H' 11 ' , analogous to the purely 
EM one in Eq. (§), 
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(49) 



and then integrated according to Eq. (Q). In a sys- 
tem of axes where the momentum transfer defines 
the z-axis and the y-axis is taken perpendicular to 
the hadronic plane defined by q and k, one has 
Wl v = H 11 + H 22 , W^'a = -2IiTLff 12 , and, in the 
"cc" treatment of current conservation, W\ v = H 00 . 
In RPWIA we use, in direct analogy with Eq. (pT|), 
the momentum-space vertex operators 



1 NC,2 — r l 7 



2m 



U = £ < e', h'\a vl , + a^57,|e, h >*< e>, h'\^\e, h >^ = ^ + - x ^ _ (fe + pf^ + ^ 



~ ^(a^s^y + iaya M „) + [/i-independent] , (47) 

where s tlv {a^ v ) are the symmetric (antisymmetric) 
parts of the leptonic tensor in Eq. (g), correspond- 
ing to the purely EM photon-exchange amplitude 
squared, |Ay| 2 |30|. Analogously, we decompose the 
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where we have not included terms corresponding to 
second-class currents An induced pseudoscalar 

axial-vector term Gpq^js is possible and will gener- 
ate a GpC7 M [. . term in the hadronic tensor. Since 
axial contributions appear only in Wy A = —2ImH 12 , 
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though, and q^ has vanishing 1 and 2 components, 
this term docs not contribute. 
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FIG. 8. Percentages of total negative-energy contribu- 
tion (1 — R ++ /R) to the parity-violating (Lav, Tav, T va ) 
responses computed in FPWIA (R ++ ) and RPWIA (R) 
with both ecl (dashes) and cc2 (dots) at |q| = 300 MeV/c. 
The A = 0.4(1) vertex functions are used in the left(right) 
panels. Notice the difference in scales. 



In the Standard Model, the electromagnetic (EM) 
and weak neutral (NC) current operators (and there- 
fore the corresponding nucleon form factors) , are writ- 
ten as different linear combinations of elementary vec- 
tor, 97^(7, and axial-vector, 97^759 (NC only), cur- 
rents. Thus, considering three flavors (q = u, d, s) 
only, and using isospin symmetry, the NC form fac- 
tors (in Sachs, G E = F 1 - tF 2 , G m = F t + F 2 , 
t = \q 2 \/4m 2 , representation) can be cast in the 
form |3S 
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V^E,M 



an /~in,v 
PV U E,M 



1 



:G 
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a — Pa u a + Pa u a 
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A 1 



(51) 



where, at tree-level in the Standard Model, 



/^ = i(l-4sin 2 



Pv — Pa — 0a — ~o 



(52) 
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FIG. 9. As in Fig. g, at |q| = 1 GeV/c. Notice the 
difference in scales. 



Provided the coefficients (3 are known (i.e., the ra- 
diative corrections small or reliably estimated) Eq. 
( |5l|) implies that independent measurements of the 
EM and NC vector current nucleon matrix elements 
allows the experimental extraction of the hitherto un- 
constrained strangeness nucleon form factors G S E M . 

Gn 
Ai 



The axial- vector form factors G^ -1 



Gj=° = G P A + G n A and G S A will also enter the parity- 
violating asymmetry A, through the Ry A response. 
Of those, G3i =0 does not contribute at tree-level, 

since, from Eq. f>|), (3 ( A ] = (3 P A + p% = °- Am7 de- 
pendence on the (pressumed small) G S A is practically 
eliminated due to the aforementioned suppression of 
R-va by oy jo- A — 10% and has to be investigated us- 
ing neutrinos |33|. Thus, only G T A =1 enters, but, as al- 
ready mentioned, the large theoretical uncertainty in 
the radiative corrections to the associated coefficient 
/?4 makes its contribution crucial. Following Refs. 
P,p3[, quantitative predictions are obtained by gen- 
eralizing the Galster parametrization to include the 
extra form factors appearing in Eq. (pT|) 



G 



T=l 



gV>{l + \ A T)-» 
G S M =h s (1 + \ s t)- 2 
G' e = Ps t{1 + X s t)- 2 

(i) 



(53) 



where, from j3 decay, at r = 0, g A — 1.26, and we use 
Xa = 3.53, A s = 4.97, i.e., equal to the Galster dipole 
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strength. Results in this section are with vanishing 
strangeness magnetic moment, /i s , and strangeness 
radius, p s §53]. 

In Figs. |8[and ^ we show the effects of the pro- 
jection on the PV responses. For brevity, we show 
results only for two choices, A = 0.4 (pv) and A = 1 
(ps), and, moreover, we do not show individual (+, — ) 
and (— ,— ) decompositions, rather only the cumula- 
tive effect of the negative energy components. The 
behavior is similar to that of the parity-conserving, 
electromagnetic responses. Specifically, the effects of 
the projection are negligible with pseudovector cou- 
pling (A = 0, not shown) for all responses and current 
prescriptions, and become more important as A — > 1; 
in all cases the effects are minimal for quasifree kine- 
matics; here, as well, the projection affects ecl more 
severely than cc2; inside the quasielastic ridge, the 
effects are more pronounced at low momentum trans- 
fer, and as with the transverse response at |q| = 300 
MeV/c in Fig. |3|, a potentially worrisome 5% effect 
is seen for the quasielastic R\\/ ecl response calcu- 
lated with pseudoscalar coupling. The only essentially 
new feature comes from the parity-violating longitudi- 
nal response, R^vi which appears extremely sensitive 
to the projection, especially towards threshold, with 
both A = 0.4 and A = 1 vertex functions at |q| = 300 
MeV/c, and, to a lesser extent, i.e., only with the 
A = 1 ones, at |q| = 1 GeV/c. This sensitivity, how- 
ever, is not a consequence of large dynamical effects 
associated with negative energy components per sr. 
Rather, it is an example of any subdominant effect 
playing an important role whenever the leading con- 
tribution vanishes (or becomes too small) , as happens 
in this case because sin 2 9w — 1/4. To see that, no- 
tice from Eq. (^l]) and the numerical values of the 
coefficients in Eq. @ w 0.046) that, although 



G E is small compared to G P E (at least for r < 1), in 
the case of the neutral current (NC) it is the neutron 
form factor, G E , that is much larger than the proton 
one, G P E . Consider then, from Eq. ( |A5| ) and Eq. ( At ), 
the ecl longitudinal response in FPWIA 
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Using Eq. ( |5l|) and Eq. (|52j), and employing the Gal- 
ster parametrization, we have 
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shell, strangeness, and T0y terms. The first tern in 
Eq. ( |55| ) arises from the electric terms in Eq. (Q), and 
the second from the magnetic terms. Because f3y is 
small, and the magnetic terms compete in Eq. (p4|), 
there is no dominant (i.e., 0(1)) term in the r.h.s. of 
Eq. (|55|). Thus, the parity violating longitudinal re- 
sponse is small in magnitude and appears with either 
positive or negative sign, depending on the kinemat- 
ics. 
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FIG. 10. Asymmetry ratios for quasifree kinematics, 
uj — uj qe . Reasults presented for two choices of electron 
scattering angles 9 and vertex functions A. Notation as in 
Fig. | 



Consider first quasifree kinematics, to = q 2 /2m. In 
this case the second term in Eq. (55) vanishes (there 
are no magnetic contributions in the "static" approxi- 
mation &). From the first term, it is clear that at low 
t values (and thus low |q|, since, for quasifree kine- 
matics, |q| = 2m \J t(1 + t ) the quasielastic PV lon- 
gitudinal response is positive and at higher r negative, 
the sign reversing at r = 0y/{lM - 5.6/3^) = 0.0278, 
or |q| ~ 317 MeV/c. The exact kinematics where 
the sign flips will depend, of course, on the small ef- 
fects not included in Eq. (p5|), as well as the projec- 
tion effects [cf. Table. |. Thus, around |q| = 320 



where we have left out small contributions from off- 



MeV/c the zeros in R\y give rise to a discontinuous 
behavior of the the ratio of R\ v in FPWIA and RP- 
WIA for quasifree kinematics, as indeed observed in 
Fig. ^. Interestingly enough, this momentum trans- 
fer value is very close to that of SAMPLE. SAMPLE, 
however, focusing on the magnetic strangeness form 
factor, is a backward angle measurement. At 8 > 50° 
the asymmetry is dominated by the T part, and, to 
a lesser extent - because of ay being small - by the 
T' part. Thus, the PV longitudinal response, and the 
projection effects associated with it, do not affect the 
SAMPLE measurement. 

Because of the behavior of R\ Vl one expects that 
projection effects will depend strongly on the electron 
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scattering angle. For quasifree kinematics, this is ev- 
ident from Fig. [k]. At |q| = 300 MeV/c, up to 4% 
effects in the forward asymmetry can be seen (this 
would be prohibitively large for strangeness studies), 
while at backward angles the asymmetry is very sta- 
ble, showing a maximum in the |q| = 300 MeV/c re- 
gion, which is, nevertheless, less than 1%. Recall that, 
in this case, the individulal R T and R\ v transverse 
responses appearing, respectively, in the numerator 
and denominator of the asymmetry, show a 0(5%) 
projection effect with ccl and A = 1 [cf. Figs. ||, || 
and even on the quasielastic peak. Due to the 
dominantly transverse character of the asymmetry at 
9 = 150°, however, these effects by-and-large cancel 
in the ratio, and at this angle the asymmetry shows 
the aforementioned < 1% variation. The extent to 
which such an effect may interfere with strangeness 
studies is examined in the next section. Fortunately, 
the projection effects become less than 1% for high 
enough momentum transfer even at forward angles, 
as, e.g., for a possible TJNAF Hall A measurement 
(|q| = f GeV/c, = 12.5°) proposed in Ref. [§. 
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FIG. 11. The PV asymmetry at |q| = 300 MeV/c 
and backward (9 = 150°) angle, using the A = 1 
Buck-Gross vertex functions. RPWIA results are shown 
with (lower) and without (upper curves) the modification 
of the on-shell form factors according to Eq. (^) . Nota- 
tion in inner graph as in Fig. |^. 



Away from the quasielastic peak, the second term 
in Eq. ( |55| ) gives a negative definite contribution, with 
a coefficient that rises with q 2 . Thus, R\ v is still pos- 
itive for low |q| (as e.g., at 300 MeV/c) and negative 
at large |q| (as e.g., at 1 GeV/c). It is always small, 
though, and, especially outside the quasielastic ridge, 
negative components constitute a significant part of 
the response, resulting in the behavior observed in 
Figs. ||and|^. At such non-quasifree kinematics, how- 
ever, R\y is quite sensitive to other aspects of nuclear 
dynamics beyond PWIA, as well, as has been shown 
in studies of both deuterium (e.g., proton- neutron in- 
terference contributions Oj) and heavier nuclei |39]. 
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FIG. 12. As in Fig. [n], at |q| = 1 GeV/c and forward 
(0 = 12.5°) angle. 



In Figs. |Tl| and [T^ we show the asymmetry in FP- 
WIA and RPWIA, as well as the corresponding ratios, 
for the SAMPLE kinematics (|q| = 300 MeV/c, back- 
ward angle) and a possible |q| = 1 GeV/c forward 
angle measurement at TJNAF. The A = 1 vertex func- 
tions are used to display the maximum effect of the 
projection. The residual relativistic effects associated 
with the relativistic treatment of the deuteron and 
the struck nucleon induce a less than 2% variation of 
the PV asymmetry within the quasielastic ridge. We 
also observe that the modification of the on-shell form 
factors to allow for a fc 2 -dependence according to Eq. 
( fi"3| ) amounts to a secondary effect (i.e., much smaller 
than the effect caused by the projection). An excep- 
tion is observed at high uo values in the |q| = 1 GeV/c 
case. This behavior is caused when deeply off-shell, 
k 2 — 0{— m 2 ), momenta, for which Eq. ( fl3| ) gen- 
erates unreasonably large modification factors, con- 
tribute to the inclusive integral in Eq. ([|). Since 
k 2 = (M d -E k ) 2 -\k\ 2 ~ 5m 2 -4my/m 2 + |k| 2 ), val- 
ues k 2 = —to 2 occur in the inclusive integral for |k| > 
1.1m. Such values are accessible for |q| > m and with 
increasing to values, since above the quasielastic peak 
(y > 0) the upper integration limit Y = y+ |q| > |q|. 

Another qualitative difference between the two 
cases, is that while at |q| = 300 MeV/c the behavior 
of ccl is quite similar to that of cc2 on the quasielastic 
peak for both FPWIA (solid lines in inner graph of 
Fig. 0) and RPWIA (dot-dashes), at |q| = 1 GeV/c 
we see (as we observed for the cross section before) 
that projection effects are larger for ccl than cc2. At 
|q| = 300 MeV/c the magnitude of the asymmetry 
is raised in going from RPWIA to FPWIA by about 
0.5 - 1.5% with either ccl or cc2 (Fig. |l|), while at 
|q| = 1 GeV/c (Fig. [l2]) there is an observable effect 
only for ccl, the magnitude of the asymmetry being 
reduced by 0.5 — 1% inside the quasielastic ridge. 
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D. Strangeness studies 

Here we discuss the extent to which the residual rel- 
ativistic effects associated with the projection affect 
the extraction of strangeness form factors from the P V 
e — d asymmetry. We are primarily interested in the 
interplay between the magnetic strangeness form fac- 
tor G S M and the NC axial isovector form factor G^ =1 
(more precisely, the radiative corrections R^ =1 to the 

relevent coefficient (3^ = 1.26[1 + R T A =1 ])- The ap- 
proved SAMPLE and E91-017 experiments will at- 
tempt to extract both quantities by combining mea- 
surements of the e — p and e — d PV backward-angle 
asymmetries. In Fig. [l3] we show "bands" of RP- 
WIA and FPWIA results for SAMPLE kinematics, 
with -R^ =1 = an d f° r three values of the magnetic 
strangeness moment fj, s entering the Ansatz in Eq. 
fl53|)) corresponding to zero strangeness (u s — 0), the 
pole model result of Jaffe (fi s = —0.3) J4(J and an un- 
likely high value (fi s = — 1). Within each /j, s "band", 
the projected (FPWIA) results show negligible vari- 
ation (0.03% on the peak) between different calcula- 
tions (i.e., ccl vs. cc2 and A = 1 vs. A = 0.4). The 
variation induced when P states and negative energy 
components are included is small, but not negligible. 
Although the ji s = 1 and fi s = bands are well 
separated, outside the quasielastic ridge the [i s = 
and (i s — —0.3 bands overlap. As hoped, though, 
they have a minimum width for quasifree kinemat- 
ics. This width depends, of course, on A, as is evident 
from Fig. ||: the A = 0.4 (long-dashes) ccl RPWIA 
curves are much closer to the overlapping upper FP- 
WIA curves compared to the A = 1 (solid) ones. On 
the quasielastic peak, the maximum variation is ob- 
served between the projected and unprojected A = 1 
ccl curves, amounting to a 0.56% effect. The asym- 
metry is linear in A([i s ) = A(0)[l — fe^/is], with 
6 P = 0.057(1) in our models. Although 6 M is approxi- 
mately tenfold the sensitivity to projection effects, it 
may not be large enough if fi s turns out to be rather 
smaller than the pole value (something not excluded 
from other calculations |if|), to the extent, of course, 
that the widths observed in Fig. |l3| give a measure of 
the theoretical nuclear physics uncertainty in describ- 
ing the e — d reaction. 

In view of the rather small sensitivity of the e — d 
asymmetry to G S E , the current consensus is to use the 
dcutcron as a means of constraining the large theoret- 
ical uncertainty in R A =1 , using primarily the proton 
for extracting G S M . Accordingly, in Fig. [l4| we fix 
the strangeness magnetic moment to the pole model 
value, and show "bands" of quasifree kinematics re- 
sults, at tree level, R A =1 = 0, as well as two other val- 
ues corresponding to the extrema in the calculation of 
Ref. §, R T A =1 = -0.34 ±0.28 [@. It is evident that 
significant constraints can be placed on R^ =1 from a 



measurement of the asymmetry and that the resid- 
ual relativistic effects discussed here should not in- 
hibit this procedure. Comparing the |q| = 300 MeV/c 
and |q| = 1 GeV/c regions, however, it appears that 
a higher q 2 measurement than that of SAMPLE (as 
for example in the E91-017 TJNAF experiment [§) 
will place tighter constraints on R A . Specifically, 
since the asymmetry is linear in Ga, we may write 
A(ns,Rl =1 ) = A( Ms ,0)[l + b A R T A =1 ]. At |q|" = 300 
MeV/c, the typical variation between models is A t h — 
0.56%, while the sensitivity to the radiative correc- 
tions b A = 0.235(4). On the other hand, at |q| = 1050 
MeV/c we find a smaller b A — 0.1259(1), but also con- 
siderably reduced model variation, A t h = 0.1%. From 
the A(n s , 0) values we find that the difference between 
the = an( l ^a =1 = —0.06 bands (upper two 

ones in Fig. ^) is 2AA th at |q| = 300 MeV/c, but 
roughly three times more effective, 7A t h, at |q| = 300 
MeV/c. 
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FIG. 13. The PV asymmetry at |q| = 300 MeV/c and 
6 — 150°. Results shown for 3 values of the strangeness 
magnetic moment /i 3 and 2 choices of Buck-Gross vertex 
functions A. For each (/i s , A) choice, results are shown with 
both ccl and cc2 for RPWIA (lower) and FPWIA (up- 
per, practically indistinguishable curves). The quasielas- 
tic peak is marked by the dotted vertical line. 



TABLE I. Momentum transfer |q cr | in MeV/c where 
Rav becomes negative (quasifree kinematics, fi s = —0.3). 
Results shown using: the static estimate in Eq. (|56|), fac- 
torized (F) ccl, and covariant (R) PWIA ccl with the 
A = 1 vertex functions. 
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In the previous section we discussed the sign- 
reversing behavior of the PV longitudinal response 
R\ v in the absence of strangeness. For quasifree kine- 
matics, where the magnetic terms cancel, only G S E 
enters, and a negative value of the strangeness radius 
modifies significantly this behavior. Using Eq. (53), 
Eq. (pa) at uo — Irtvr now becomes 
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The resulting quadratic equation gives an estimate 
of the momentum transfer value |q CT | where R\ v re- 
verses sign as a function of p s . As seen from Table 
[jj this "static" estimate is reasonably accurate. For 
positive strangeness radius |q cr | depends mildly on 
p 8 , but for negative values |q cr | becomes significantly 
larger and for p s < —1.4 (the pole model value is 
—2.1 ± 1.0) the response is always positive. 

Thus, given the large variation in the model predic- 
tions of p s pC| ], a determination of R\y m different |q| 
values might be helpful in constraining p s . Although 
quite challenging, a separation of the L, T and T' 
parts of the PV asymmetry is not beyond current ex- 
perimental capabilities, and will be performed for the 
e - p PV asymmetry as part of the E91-017 TJNAF 
experiment [p3|. 
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FIG. 14. The PV asymmetry [in units of q 2 ] as func- 
tion of momentum tranfer |q| for quasifree kinematics and 
9 = 150°. Results shown for 3 values of fij^ (see text). 
Solid lines: ccl, dashes: cc2. A magnification of the 1 
GeV/c region is shown in the inner graph. 



V. SUMMARY 

In this work parity-violating, electron-deuteron 
scattering in the quasifree region has been studied in 
relativistic PWIA. The emphasis has been to quantify 
the role of residual relativistic effects associated with 
the non-relativistic treatment of the bound state in 



the semirelativistic, factorized PWIA, assuming that 
the latter correctly incorporates relativistic effects as- 
sociated with the large momentum transfer. The main 
results of this investigation can be summarized as fol- 
lows: 

• In the covariant calculation (RPWIA), factor- 
ization is spoiled by off-shell contributions as- 
sociated with the struck nucleon. A weak fac- 

(56) torization still holds, with the momentum dis- 
tribution and the single-nucleon tensor being 
now density matrices in spin and energy-sign 
spaces of the struck nucleon. The familiar FP- 
WIA using de Forest's single-nucleon tensor cor- 
responds to a projection of the RPWIA to the 
positive energy sector. That is, the -ff++ single- 
nucleon tensor is the FPWIA single-nucleon ten- 
sor, and the p++ momentum distribution is 
the nonrelativistic deuteron momentum distri- 
bution. No additional operations (such as ne- 
glect of P states that superficially appear in 
p++, or nonrelativistic expansion) is required. 
Thus, this truncation defines the residual rel- 
ativistic effects we are interested in, and shows 
the close connection between relativistic and off- 
shell effects. Similar considerations apply to 
quasielastic scattering on any nucleus and not 
just deuterium. 

• Since, from the above, a covariant treatment 
of the e — d scattering requires inclusion of off- 
shell effects, it is not consistent to incorporate 
such effects in the d-n-p vertex but ignore them 
in the jNN and Z°NN vertices. A measure 
of the thusly induced model dependence is ob- 
tained by using various on-shcll-cquivalcnt pre- 
scriptions for these vertices, and several model 
calculations for the deuteron vertex functions. 

• The residual relativistic effects depend strongly 
on the nNN coupling used for calculating the 
deuteron vertex functions (minimal for pseu- 
dovector and maximal for pseudoscalar cou- 
pling). They also depend on the off-shell pre- 
scription for the nucleon vertex. Although the 
difference between ccl and cc2 prescriptions is 
negligible in FPWIA, it is not such in RP- 
WIA, where projection effects appear more pro- 
nounced with the ccl prescription. In that re- 
spect, the cc2 presription seems preferable. 

• Effects associated with allowing for an off-shell 
dependence of the nucleon form factors are 
generically smaller than projection effects. The 
arbitrariness of the Ansatz used, however, does 
not allow definitive conclusions to be drawn. 

• With the exception of the parity-violating lon- 
gitudinal response, R\ v -, inside the quasielastic 
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ridge the residual relativistic effects are below 5- 
10% with ccl, and a few % with cc2, when the 
pseudoscalar coupling (which maximizes such 
effects) is used. On the quasielastic peak, ef- 
fects are typically below 2%. A notable excep- 
tion is the |q| < 400 MeV/c region where the 
transverse responses, R T and R^ v , can have 
> 4% effects. The R h AV response can show quite 
important effects, especially in the |q| = 300 
MeV/c region, which reflect the lack of a dom- 
inant term in this response, thereby exposing 
otherwise non-leading contributions. 

The magnitude of these corrections for quasifree 
kinematics is small for them to play a role 
in the cross section. In the parity-violating 
asymmetry, where 1% effects are important for 
strangeness studies, the effects mostly cancel in 
the ratio. The R\ v response is sensitive to these 
effects and, thus, effects > 2% are seen in the 
forward angle asymmetry for |q| < 800 MeV/c. 
For backward scattering, where the deuteron 
measurements will be used for disentangling the 
correlation between magnetic strangeness and 
axial-vector isoscalar radiative corrections, the 
PV asymmetry is sensitive to less than 1% to 
residual relativistic effects, the maximum sensi- 
tivity oobserved in the |q| = 300 MeV/c region. 
Thus, the extraction of G^ =1 at the E91-017 
TJNAF kinematics will be less sensitive to such 
effects than at the SAMPLE kinematics. 
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APPENDIX A: 

Here we present some details of the calculation of 
the hadronic tensor in RPWIA. After defining pro- 
jected vectors 



V = r" - d 



W = 



Ml 
M3 ' 



(Al) 



a straightforward calculation gives for the scalar term 
P in Eq. (|) 



R 2 A R 

/? = 3mA 2 + — V 2 -2— (V-n) 
m m 



AF 
BG 



3m- 3- + 2^ '- 



-V 1 



1 + 



m 
(k-n) 



m 

AG + BF 

m 



[{k-V) + (n-V)] 



G 2 
4m 3 
F 2 



V 2 [( 



2 |7™2 , h 2\ 



2(k-n)\ 



+ U(k-N) - 6(k-n) + 3(m 2 + k 2 )] 

4m L J 

FG 

+ ^ [-2(m 2 + k 2 )(V-n)~4m 2 (k-V)] , (A2) 
and for the vector term 



V = A 2 [3n M - 2N"] - 



B 2 V 2 



2ABV^ 



AF [3^ - 2N» - 3fc^] 

^rV 2 [Jb" + n"] 
m z 

AG + BF \ m 2 V » + (n-VW] 
m 

AG-BF 



[(k-n)V^ L - (k-V)n^] 



Tfl 

+ [(fc 2 - m 2 )n» - 2 [(k-n) + m 2 ] F] 

F 2 

+ [[6(k-n) - A(k-N) - 6 m k 



FG 

4m 3 



+ (k 2 - m 2 ) [2N V - 3n u ] 
[4m(P-n)fc ) ' + 4m(tI))F 
+2m(m 2 - k 2 )V] . 



(A3) 



Using the f£ vertices of Eq. @ and Eq. (jlj), the 
cc2 interference hadronic tensor [cf. Eq. (|4^)] reads 

Z H£ 2 = F l F l [S» P V + Sy + {mT - (S-p))] 

+ ^r 1 W ( T M - m ( s -i)) + 1^ ( mSV - t p v )] 

2m 

+ ^ [<f (T(p.q) - m(S-q)) + q v (mS» - Tp»)} 

+ [9^ W{S-p) + mTq 2 - 2(S-q)(p-q)] 

-q»q v [(S-p) + mT] + (p-q) (S»q v + S v q») 
+ {S-q) (p^q" +?Y) - q 2 {S»p v + S"/)] 

+ tG A F 1 e^ ab S aPb + i^^e^ ab q b (mS a + Tp a ) , (A4) 
2m 

with Z = 3(fc 2 — m 2 ) 2 /4. Analogously, for ccl 
Z = G m Gm [W + S>" + gT (mT - (S-p))} 
-(k+pY(mS u +Tp v ) 

'-(k + p) v (mS" + Tp") 



F2GM . 

2m 
GmF 2 fT 

2m 
F2F2 



(k+pY(k+pY (mT+(S-p)) 



Am 2 

+ i G A G M e^ ab S aPb 



(A5) 



1 « 



The ot her tensors us ed i n this work can be mapped 
to Eq. and Eq. (|A5|) 



• for the electromagnetic hadronic tensors, Pi 2 — ► 
Pi,2, Ga - 0. 

for the corresponding FPWIA tensors of Eq. 
I5|), 



(A6) 



The ccl [Eq. (|Aj|)] and cc2 [Eq. (A4)] forms agree in 
the "on-shell" limit, consisting of (I) Eq. (A6), and (I) 
letting — > <p, where = p^ — fc M . For example 



under Eq. (A6) the P1P2 terms in Eq. (A4) and Eq. 
jAq) become 

cc2^~[^(g- g )-^r] 

ccl -> - [2.g^" y (m 2 - k ■ p) - (k - p)*(k - p) v ] 

= [9^{q-q)~q»<f] , (A7) 



where we have set c = f + f^/12ir 3 AId(k 2 — m 2 ) 2 , and 
used P4 7 ° = (p+ $/2, P fe "P° = (k-V + n-V)/2, and 
£| = (to 2 + fc-nj/2. We finally obtain 

(p+-,rs'H+--rs + P-+;rsW W + . rs ) 

spins 

U 



/+/- [m 2 



f_ 

h ' /+ 
(fc-n)] 



//.;/ 



127r 3 M d (fc 2 -TO 2 ) 2 



Tr 



2to 



2m 



(B3) 



Let us compare Eq. (Be) with Eq. (|17|). The first two 
"mixed" terms in Eq (B3) correspond to the (+, — ) 
contribution to the first term in the RHS of Eq. (p7[), 
while the last term, which spoils factorization, derives 
from the second term in the RHS of Eq. ([l7|), as an- 
ticipated in Eq. (Bfl). 



which coincide under q^ — ► q^ 



APPENDIX B: 



Here we evaluate the (+,— ) and (— ,+) contribu- 
tions to Eq. (^3|), using the relevant spectral function 
[cf . Eq. (]2q)] . Notice that we have chosen to define £ M 
in terms of the positive energy spinors u(k, s). Thus, 
k + -£_ = but = n-£ = 2E k £° ^ 0. The spin sums 
are evaluated_using J2 S £L = and J2 S J2 r = 

2 (~<7 Ml/ + fc^Tc^/m 2 ). Since p_| is linear in £ M and 

because of the first of these properties, there is no 
contribution in Eq. (|2^) from the <5 s>r part of the pro- 
jector in Eq. (p6[). Thus, effectively, 



x 8 +x r - + x 8 -x+ 



75 ( 757 + 757 



^+TO 
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(0) ~ 



S5,r 



4m 4m 
= (E k r - (p+m)r°) /2m , 

Therefore, from Eqs. (E 
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(Bl) 
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